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Abstract 

\ The electromagnetic pion form factor, F n (q 2 ), is calculated for spacelike-c/ 2 in 

impulse approximation using a confining quark propagator, S, and a dressed quark- 
photon vertex, r„, obtained from realistic, nonperturbative Dyson-Schwinger equa- 
tion studies. Good agreement with the available data is obtained for F n (q 2 ) and 
other pion observables, including the decay 7r° — ► 77. This calculation suggests 
that soft, nonperturbative contributions dominate F 7T (q 2 ) at presently accessible q 2 . 
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1 Introduction 



As a bound state of a light quark and antiquark, the pion is an ideal system for 
exploring the application of different approaches to the study of bound state 
structure in QCD, which is intrinsically nonperturbative. Such studies are 
constrained by the Goldstone boson nature of the pion. The internal structure 
of the pion affects its observable properties. The pion electromagnetic form 
factor, F v (q 2 ), is one observable that is sensitive to this internal structure and 
it has been much studied (JTJ 

Perturbative QCD has been employed in order to estimate the behaviour of 
F n (q 2 ) at large spacelike-g 2 . These analyses rely on the separation of the 
amplitude into a product of soft and hard contributions using a factorisa- 
tion Ansatz [ffl], however, the applicability of this approach to exclusive pro- 
cesses is uncertain. In this factorisation approach F n (q 2 ) is the product of a 
soft contribution that depends on the bound-state Bethe-Salpeter amplitude, 
which provides only the overall normalisation, and a hard contribution that 
is independent of the bound-state Bethe-Salpeter amplitude and is taken to 
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be given by the Born amplitude for a collinear quark- ant iquark pair, each 
massless, to scatter coherently from a virtual photon. This analysis yields 
q 2 F n (q 2 ) = 16 a(q 2 ) f 2 as spacelike-g 2 — > oo, where f n ps 92 MeV and a(q 2 ) is 
the running coupling constant in QCD. 

It is not clear whether presently accessible values of q 2 are large enough to 
test predictions based on perturbative analyses in QCD. It is argued in Ref. |J 
that they are not; i.e., that the factorisation Ansatz is invalid at presently ac- 
cessible values of q 2 and hence that the q 2 dependence of the quark momentum 
distribution in the pion provides an important contribution to F n (q 2 ). This 
conclusion is supported by Ref. || and by the fact that a good fit to the 
experimental data, over the entire range of available q 2 , is possible using the 
light-front formulation of a relativistic constituent quark model |IJ , which has 
no obvious connection with perturbative QCD. 

Herein the impulse approximation to F^(g 2 ), illustrated in Fig. [I], is calcu- 
lated for spacelike-g 2 as a phenomenological application of the nonperturba- 
tive Dyson- Schwinger equation [DSE] approach to QCD, which is reviewed in 
Ref. 0. The primary elements of this calculation are: 1) The dressed quark 
propagator, S(p), which is confining in the sense that it has no singularities 
that can lead to free-quark production thresholds in Fig. [IJ i.e., there is no 
quark mass-shell; 2) The pion Bethe-Salpeter amplitude, T n (p,P), which is 
regular for spacelike values of p, the relative q — q momentum. In the chiral 
limit, m — > 0, T w (p,P) is completely determined by the dressed quark prop- 
agator, which is a manifestation of Goldstone's theorem in DSE approach ||; 
and 3) The dressed quark-photon vertex, T fl (pi,p 2 ), which is regular in the 
spacelike region; i.e., away from resonances such as the p-meson, and follows 
from extensive QED studies [fHTTH. These properties, which together are suffi- 
cient to ensure confinement, entail that this impulse approximation calculation 
is free of both endpoint and pinch singularities, which arise in perturbative 
analyses. 



In phenomenological applications of the DSE approach the model dependence 
is restricted to spacelike-g 2 < 2 GeV 2 , and is realised in a modelling of the 
form of the quark-quark interaction in the infrared. This not only incorporates 
information obtained about, for example, the gluon condensate in the QCD 
sum rules approach [12| but also extends it. The calculation of experimental 
observables in this approach therefore allows one to place constraints on the 
qualitative and quantitative features of the effective quark-quark interaction 
at small spacelike-g 2 in QCD and to infer the q 2 scale where perturbative, 
model- independent, effects begin to dominate. 



In Sec. 2 the impulse approximation and its primary elements [S, T n and V^] 
are discussed in detail. The width T 

71*0- — ^'V'y IS calculated in impulse approxima- 
tion in Sec. 3 and shown to be independent of the details of S, T n and in 
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Fig. 1. A pictorial representation of the amplitude identified with the tt-tt-A^ 
vertex in impulse approximation: pions - thick, straight external lines; (Tr\qq) 
Bethe-Salpeter amplitudes - filled circles at the 7r-legs; photon - wiggly line; dressed 
quark-photon vertex, which satisfies the Ward-Takahashi Identity, Eq. ( p4| ) - shaded 
circle at the the photon-leg; dressed quarks - thin internal lines. 

the chiral limit. This illustrates the manner in which anomalies are realised in 
the present framework. The calculation of F n (q 2 ) for spacelike-g 2 , and other 
pion observables, is described in Sec. 4 and the results compared with experi- 
ment. The behaviour of F n (q 2 ) at large spacelike-g 2 in impulse approximation 
is determined analytically in Sec. 5: F^q 2 ) oc l/q 4 , up to ln[g 2 ]-corrections. 
This result is verified numerically and found to become dominant only for 
spacelike-g 2 > 10 GeV 2 , which is presently inaccessible experimentally. This 
asymptotic form is a consequence of the realisation of confinement explored 
herein. Its validity or otherwise is not an essential consequence of the DSE 
framework. The results are summarised and conclusions presented in Sec. 6. 



2 Impulse Approximation 

Herein all calculations are carried out in Euclidean space, with 7 M hermitian 
and metric 5^ u = diag(l, 1, 1, 1). 

One may define the impulse approximation to the connected n-n-A^ vertex in 
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QCD as, with m u = 



2N f r— 

A M (P + q, -P) = / (0 ti D [T^k; P + q)x 

S(k ++ )ir„(k ++ , k. + )S(k_ + )T n (k - \q- -P)S(k^ 

where q is the photon momentum and P is the initial momentum of the pion. 
Here the trace over colour and flavour indices has been evaluated leaving only 
the trace over Dirac indices and 

Kp = k + § q + f P . (2) 

In Eq. (HI): T fl (pi,p 2 ) denotes the dressed quark-photon vertex; T n (p;P) the 
pion Bethe-Salpeter amplitude, with p the relative momentum and P the 
centre-of-mass momentum; and S(p) the dressed quark propagator. The Bethe- 
Salpeter wave-function is 

xAp; p) = sip + \p) r ff (p ; p) s( p - \p) . (3) 

The impulse approximation, Eq. ([I]), is illustrated in Fig. [I] and can be derived 
as an application of the formalism described in Ref. Its definition is only 
complete when the functions S, T n and have been fully specified, which is 



the subject of Sees. |2T2| and respectively. 
2.1 Quark Propagator 



The dressed quark propagator in Eq. ([[]) can be obtained by solving the fol- 
lowing Dyson- Schwinger equation [DSE]: 

S- 1 (p)=z 1 -p + m+lg 2 J ^^S(k)n(k,p)D^(p-k), (4) 

where m = (m u + md)/2 is the current-quark mass. Here D^ u {k) is the dressed 
gluon propagator and r^(pi,p 2 ) is the dressed quark-gluon vertex; each of 
which satisfies its own DSE. The general form of the solution of Eq. is 

S(p) = -i-y -pa v (p 2 ) + cr s (p 2 ), (5) 

which can also be written as 

S ^ = ~ A( 2\ l 1 U ( 2\ ■ ( 6 ) 

• pA{p z ) + m + B{p z ) 

This equation has been much studied and the general properties of as and cry 
in QCD are well known QJ. 
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As a simple example, in Ref. [M| Eq. ([|) was solved with 



9 2 D fxu (k) = ( ^ - ^ ) 8<K 4 D5 4 (k) and r M (p,p) = -ifyS^Cp) , (7) 



where D is a mass-scale parameter. This Ansatz for the dressed gluon prop- 
agator models the infrared behaviour of the quark-quark interaction in QCD 
via an integrable infrared singularity, as suggested by Refs. |15|-|l7| , and and is 



sufficient to ensure confinement, in the sense described below. The Ansatz for 
the dressed-quark gluon vertex is the result of extensive analysis of its general 
form HU. The solution of Eq. (|) in Ref. is 



<rs(v) = 
C ^2 



exp(-2y 2 )J 1 (Arhy) + — / d££ #i(m£) Hvt) e*P 

v J 



2my y 

u 

with m = m/y2D, y 2 = p 2 /(2D) and where J\ and K\ are Bessel functions, 
and 

°v(y) = -(crs(y) + ^-^-cTs(y)) , (9) 

m \ Ay ay J 

with a s (y 2 ) = V2Da s (k 2 ) and a v (y 2 ) = 2 D oy{k 2 ). In Eq. (|), C is a 
parameter associated with dynamical chiral symmetry breaking and it is not 
determined by Eq. @) with Eq. (0), while the integral, which is associated 
with explicit chiral symmetry breaking, cannot be evaluated in terms of known 
functions. 

At large spacelike-p 2 one finds from Eqs. (^) and (|9]) that 



. 2 . m m 3 , . 2 . 1 D + m 2 
°s(p 2 --- — + ••• and cr V (p 2 )~ — — + .... (10) 

pZ pi pZ pi 

At large spacelike-p 2 in QCD one has at leading order 



/;/ 



°s{p 2 ) w — — d (11) 



p2 



with m a renormalisation point invariant and d = 12/ [33 — 2Nf]\ Nf is the 
number of quark flavours. One therefore sees that, neglecting ln[p 2 } terms, the 
model defined by Eqs. (|]) and (0) incorporates asymptotic freedom. 

Another feature of this model is that both as and ay are entire functions in 
the complex-p 2 plane with an essential singularity. As a consequence the quark 
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propagator does not have a Lehmann representation and can be interpreted 
as describing a confined particle. This is because, when used in Eq. ([!]), for ex- 
ample, this property ensures the absence of free-quark production thresholds, 
under the reasonable assumptions that T n is regular for spacelike-p 2 and 
is regular for spacelike-g 2 . It follows from this that Eq. ([I]) is free of endpoint 
and pinch singularities. 

This is a particular, sufficient manner in which to realise the requirement that 
Fig. [l] have no free-quark production thresholds, which is the definition of 
confinement explored herein. There are other, more complicated, means of 
realising this definition but the effect is the same [0] . Some of the phenomeno- 
logical implications of a model with a simple realisation of this confinement 
mechanism have been discussed in Ref. HIT 



The solution described by Eqs. (|8|) and (||) has a defect. To see this one sets 
fh = in Eq. (||]), which yields 

a s (y) = Ce- 2 ^ . (12) 
This poorly represents as for sufficiently large y since, in the absence of a bare 



mass and when chiral symmetry is dynamically broken, it is known [19] that 



4ir 2 d K 

*g(y)Uo-> — y2(lny2) i- rf ( 13 ) 

with k = — (ln[n 2 /AQ CD ])~ d (qq)^2, a renormalisation point invariant. This 
defect results from the fact that, although the form of D^k) in Eq. (|7|) 
generates confinement, it underestimates the strength of the coupling in QCD 
for sufficiently large k 2 . In the numerical studies of Eq. that have used a 
better approximation to D^ u {k) [^Oj there is no such defect. 



2.1.1 Approximate, algebraic quark propagator 

Herein, to avoid the need for a numerical solution of Eq. (^), Eqs. and 
d^) are simply modified so as to restore the missing strength at intermediate- 
x (= y 2 ) and thereby provide a better approximation to the realistic numerical 
solutions p0| , while retaining the confining characteristics present in the model 
example described in association with Eqs. (f|) and (|). 

The following approximating algebraic forms are used: 



1-e 

b\X 



= Ce~ 2x + 

-bix J _ e ~b 3 x 



(14) 



b 3 x 



bo + k 



1-e 
Ax 



-Ax 1 



+ 



rn 



x + m 2 



— e 



-2 (x+fh 2 y 
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2(x + m 2 ) - 1 + e - 2 ^+™ 2 ) _ _ 2x 
°v{x) = - ; _ 2 , 2 mCe 2x , 15 

which are entire functions, as in the model described above, and allow the 
propagator to be consistent with realistic numerical solutions of Eq. @. 

When 60 = = 62 , Eqs. ([14]) and (|1^) provide an excellent approximation to 
Eqs. (U) and @ ||21|| , while for nonzero values of these parameters it is clear 
that the behaviour given in Eq. ( Jl3| ) is recovered, up to ln[p 2 ]-corrections. 
These model forms are also entire functions in the complex p 2 plane with an 
essential singularity. 

The expressions in Eqs. ([L4]) and (|15D provide a six parameter model of the 
quark propagator in QCD: C, m, b , . . . , 63. [A(= 1CT 4 ) is introduced simply 
to decouple 62 from the quark condensate, as will be shown below.] These pa- 
rameters can be fitted to experimental observables and Eq. used to place 
constraints on D tiu (k), the effective gluon propagator. In S one therefore has 
an implicit parametrisation of D^ik) and hence a connection between exper- 
imental observables and the nature of the effective quark-quark interaction 
in the infrared. A phenomenologically successful application of this procedure 
may then make possible the use of precise experimental data as a probe of the 
effective quark-quark interaction in the infrared. 



2.2 Pion Bethe-Salpeter Amplitude 



The Bethe-Salpeter amplitude in Eq. (JT|) is the solution of the homogeneous 
Bethe-Salpeter equation [BSE]: 

K S (P;P) = J $jtK™*»(k,p;P) (S(k + lP)T 7T (k;P)S(k-±P)) tU (16) 

where P is the centre-of-mass momentum of the bound state, p is the relative 
momentum between the quarks in the bound state and the superscripts are 
associated with the Dirac structure of the amplitude. In the isospin symmetric 
case, m u = rrid, K rs ' tu (p, k; P) oc If, the identity matrix in flavour-space. 
Further, since T n and S(p) are oc Ic, the identity matrix in colour-space, then 
K rs ' tu (p,k;P) oc I c also. 

The generafoed-ladder approximation is defined by the choice 

K^ tu (p } k;P) = l g 2 D, v (p - k) { lv ) us (17) 

in Eq. fllED, with D^ u (p — k) the dressed gluon propagator and S(p) the dressed 
quark propagator obtained from Eq. (f|). This equation has been much stud- 
ied [^] . In using dressed quark and gluon propagators, the generalised- ladder 
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approximation is a significant advance over the ladder approximation familiar 
in QED where perturbative fermion and photon propagators are used. The 
amplitude obtained as a solution in generalised-ladder approximation defines 
the nonperturbative "dressed-quark core" of the meson and provides the dom- 
inant contribution to physical observables, as will be seen herein. 

The most general form of allowed by Lorentz covariance, which is odd 
under parity transformations, is 



T n (p;P)= il5 ^E(p;P)+ (18) 

i n . pp. PF(p; P) + z 7 • PG(p; P) + [7 -p, 7 • P] H(p; P)\ 

where, since 7r° is even under charge-conjugation, E, F, G and H are even 
functions of (p ■ P). The many studies of Eq. ( |TE| ) using Eq. ( |TTD |2"2"| ] suggest 



that the dominant amplitude in Eq. (^) is E(p; P), with the other functions 
providing ~ 10% contribution to physical observables; i.e., that it is a good 
approximation to write 

T n (p,P)=t l5 E(p;P) . (19) 

In generalised-ladder approximation K is independent of P and the standard 
normalisation condition for T n reduces to the requirement that, for P 2 = —m 2 : 



x tro 



r.(fc; P)^^-Y n {k- -P)S(k -) + 



T 7T (k;P)S(k 0+ )T n (k;-P) 



dP, 
dS(k .) 
dP„ 



(20) 



2.2.1 Pion as Goldstone mode and quark- antiquark bound state 

In the chiral limit; i.e., when the current quark mass, m, is zero, the pseu- 
doscalar generalised-ladder approximation BSE and quark rainbow-DSE [which 
has = 7^ in Eq. @] are identical |§. This entails that there is a massless 
excitation in the pseudoscalar channel with 

E(p;P 2 = 0) = ^B m=0 (p 2 ) , (21) 

where B(p 2 ) is given in Eq. (Q) and is the calculated normalisation constant. 
This is the manner in which Goldstone's theorem is realised in the Dyson- 
Schwinger equation framework. In this case, Eqs. (|14|) and (|15|) completely 
determine T n . 
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The dichotomy of the pion as both a Goldstone boson and a quark-antiquark 
bound state is thus easily understood in the DSE approach. One has dy- 
namical chiral symmetry breaking when, with m = 0, Eq. (|J) yields a solution 
B m=0 (p 2 ) ^ 0; i.e., when a momentum-dependent quark mass is generated dy- 
namically by the interaction of the quark with its own gluon field. For m = 
one also finds that the Bethe-Salpeter equation in the pseudoscalar channel 
reduces to the quark DSE as P 2 — ► 0, where P M is the total-momentum of the 
bound state. It therefore follows, without fine tuning, that if the quark-quark 
interaction is strong enough to support dynamical chiral symmetry breaking 
then one has a massless, pseudoscalar bound state of a strongly dressed quark 
and antiquark whose Bethe-Salpeter amplitude is the quark mass function, 
Eq. (pi]) ; i.e., one has a zero mass bound state of a quark and antiquark, each 
of which has an effective mass of ~ 200-400 MeV. (The actual value depends 
on the definition and the strength of the interaction.) 

For m^0, the pion Bethe-Salpeter amplitude must still vanish as the relative 
momentum p 2 — > oo A first, simple approximation in this case is 



E(p; P 2 = -ml) « 1 B m ^{p 2 ) , (22) 

JTT 

which, for small current-quark mass, is very good both pointwise and in terms 



of the values obtained for physical observables f2l) |. One notes that using 
Eqs. (HD and (pD in Eq. (f|) entails 

r w (p,P)oc^, (23) 

which, up to the ln[p 2 ]-corrections associated with the anomalous dimension, 
reproduces the ultraviolet behaviour of the Bethe-Salpeter amplitude given by 
QCD [M. 



2.3 Quark-photon Vertex 



The quark-photon vertex, T^(pi,p 2 ), satisfies a DSE that describes both strong 
and electromagnetic dressing of the interaction. Solving this equation is a 
difficult problem that has recently begun to be addressed |26[]. Much progress 
has been made in constraining the form of T fl (p 1 ,p 2 ) and developing a realistic 



Ansatz [1C,11 



The bare vertex: T^(pi,p 2 ) = 7 M , is inadequate when the fermion propagator 
has momentum dependent dressing because it violates the Ward-Takahashi 
identity: 



(Pi - P2)^ln ^ S (pi) - S (pa) 



(24) 
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and hence leads to an electromagnetic current for the pion that is not con- 
served. 



An Ansatz [[H]] fulfilling the criteria [11 that it a) satisfies the Ward-Takahashi 
identity; b) is free of kinematic singularities; c) reduces to the bare vertex in 
the free field limit as prescribed by perturbation theory; and d) has the same 
transformation properties as the bare vertex under charge conjugation and 
Lorentz transformations, is 



(25) 



where 



P 



k 2 



\A{p 2 ) + A{k 2 )} 



A{p 2 ) - A(k 2 



[7 • p + 7 • k] 



(26) 



— 1 



B{p 2 ) - B(k 2 



and (p—k)^ Tj(p, k) = with rj(p, p) = 0. r^ c in Eq. fl2l>D is completely deter- 
mined by the dressed quark propagator but must be determined otherwise. 
Gauge covariance of the fermion propagator and multiplicative renormalisabil- 



to place constraints on Tj^(p,k). 
1 and B = constant, TT, = 0. 



ity of the fermion DSE can be used 
Using the bare quark propagator, which has A 

The definition of the impulse approximation is complete with 

r»( P ,k)^T* c (p,k) 



(27) 



in Eq. d). 



2.3.1 Vector meson dominance 

The electromagnetic pion form factor, F n (q 2 ), is an analytic function but for 
a cut extending from q 2 = —00 to —4m 2 ; i.e., on the timelike real axis. It is 
real for q 2 > —4m 2 . Generalised vector meson dominance is the observation 
that the pion form factor is completely determined by the spectral density 
associated with the J PC = 1 pion-photon vertex, which is the discontinuity 
across this cut. 

Models of vector meson dominance consist in a phenomenological Ansatz for 
this spectral density. A common model is a single, p-meson pole representa- 
tion. Since the vertex in Equation (p7|) is free of kinematic singularities, by 
construction in accordance with criterion b) above, vector meson pole contri- 
butions of this type are explicitly excluded in impulse approximation. 
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Explicit vector- meson-photon mixing contributions to F n (q 2 ) enter through 
the dressed quark-photon vertex and only appear in T^(p,k) because vector 
mesons are resonances, which are only uniquely defined on shell where they 
are transverse. They appear via resonance contributions in the timelike region 
to the hadronic component of the photon vacuum polarisation. 

For q 2 = the differential form of the Ward identity, T^(p,p) = —id^S(p), 
completely determines the quark-photon vertex (r^(p,p) = 0). This entails 
that the pion charge radius is primarily determined by the dressed quark 
propagator and is not very sensitive to T^(p,k). In the sense of generalised 
vector meson dominance described above, this admits an interpretation that 
the charge radius is primarily determined by the non-resonant part of the full 
spectral density near threshold. To go further and identify this as the "tail 
of the p-meson resonance" involves additional model assumptions since the 
vector mesons are resonances, whose definition off-shell is arbitrary. 



This point is also discussed in Ref. [27| 



2.4 Current Conservation 



Using Eq. (^7]) and the identities: 

S(-k) T = Ct S{k) C- r£(-fc; p) = Ct T^k; p) c- 
Tli-k; p) = ^ T 7T (k; p) C; T T (-k, -p) = -C* T,(p, k) C, 



(2* 



where C = 7274 is the charge conjugation matrix, one finds easily that the 
7r-current is conserved: 

g M A M (P + q,-P)=Q. (29) 
For elastic scattering, with [2q ■ P + q 2 = 0] , one can therefore write 

A M (P + q, -P) = (2P + q)^ F„(q\ P 2 ) ■ (30) 
One obtains similarly that 



A M (P,-P) = 2P^F W (Q,P 2 



2N r 



d 4 k 



tr D 



r^;P)^-x 



BP, 



T^k; -P)S(k _) + T T (k; P)S(k 0+ )T 7T (k; -P) 



dS(k _ 
dP tl 



(31) 



Comparing this with Eq. (B^) it is clear that in impulse approximation F n (0, P 2 ) 
1 only if the Bethe-Salpeter kernel is independent of P; i.e., impulse ap- 
proximation combined with a P-independent Bethe-Salpeter kernel provides 
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a consistent approximation scheme. In this case one has, in the chiral limit 
P 2 = -ml = 0: 



j - 



Nr. 



ds s B(s) 



(32) 









ii I l \2 




s 




-s 2 




) 



with s = p 2 



This shows that the impulse approximation to the form factor, Eq. ([j]), is 
regular in the chiral limit. In fact, the calculated results are only weakly de- 
pendent on m„ 



The study of Ref. §27 



indicates that, at m n = 0.14 GeV, 
Eq. (JX|) provides the dominant contribution to F n (q 2 ) for spacelike-g 2 and that 
pion-loops are unimportant. 



3 Neutral pion decay 



In Euclidean space the matrix element for the decay n — > 77 can be written 

ex 

M{h,k 2 ) = -2i-^e tiUpr7 e fl (k 1 )e l/ (k 2 )k lp k2 r7 G(k 1 - k 2 ) , (33) 

where ki are the photon momenta and e(fcj) are their polarisation vectors. 
Here, the ir° momentum is P = (ki + k 2 ) and P 2 = 2 k\ ■ k 2 . 

Using Eq. ([33|) one finds easily that 



,3 



2 



r.o^ 77 = ^ ) G(-ml) 2 . (34) 

16tt 

Experimentally one has r^o^^ = (7.74 ± 0.56) eV, which corresponds to 

^o 77 = G(-ml) = 0.500 ± 0.018 , (35) 
using m^o = 135 MeV and f w = 92.4 MeV. 
In impulse approximation one has 



i—. - < avoa hp k 2a G(h ■ k 2 ) = J d 4 k ti D 



S(k) (36) 



a, 

iT^k, k - h) S(k - h) T^k + k 2 /2- P) S(k + k 2 ) iT u {k + k 2 , k) 
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In the chiral limit, P 2 = 0, and using Eq. ([19]), one easily obtains 



3°o 77 = G(0) = J dss [U E(s; 0)] A a v x (37) 
o 

[A [a v a s + s (a' v a s - a v a' s )) + sa v (A' a s - B'o v ) J . 
Now, denning 

C(s) = ^ = ^ (38) 



one obtains a dramatic simplification and Eq. (|37|) becomes 

o° - 7dc UE{3]0) C ' {S) (39) 

It therefore follows from the manner in which dynamical chiral symmetry 
breaking is realised in the Dyson- Schwinger equation framework; i.e., from 
Eq. (pip , that in the chiral limit 
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*" = /'" (T+cF = 5' (40) 

v y 

since C(s = 0) = oo and C(s = oo) = 0. Hence, the experimental value is 
reproduced independent of the details of S(p). 

This illustrates the manner in which the Abelian anomaly is incorporated in 
the Dyson- Schwinger equation framework. Similar results are obtained for the 



Wess-Zumino five-pseudoscalar term p|| and 'jtitttt interaction [29|. 



In order to obtain the result in Eq. fllOD it is essential that, in addition to 
Eq. the photon-quark vertex satisfy the Ward identity. This is not sur- 
prising. However, the fact that one must dress all of the elements in a calcu- 
lation consistently is often overlooked. 

The subtle cancellations that are required to obtain this result also make it 
clear that it cannot be obtained in model calculations where an arbitrary 
cutoff function (or "form-factor") is introduced into each integral. The fact 
that E(p; P 2 = 0) is the pion Bethe-Salpeter amplitude and f n E(p; P 2 = 
0) = B(p 2 ) in the chiral limit, is crucial. 

A good approximation to g^yy for nonzero pion mass is given by Eq. (|36|) 
with Eq. <m). 
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4 Calculated Spacelike Form Factor 

The impulse approximation to F n (q 2 ) is denned by Eq. (P with the dressed 
quark propagator obtained from Eq. (|]), the Bethe-Salpeter amplitude ob- 
tained using Eqs. flUf), (|17D and the dressed quark-photon vertex in Eq. (p7|) . 

If, at large spacelike-g 2 , each of the elements in Eq. ([[]) behaves as prescribed 
by the renormalisation group in QCD then this amplitude reduces to that 
studied in perturbative analyses of the elastic form factor. [To see that the 
"hard-scattering" contribution is included one need only once-iterate the ir- 
q-q vertices using Eq. fllBf).] It is therefore plausible that a model such as 
the one described by Eqs. (|14"D , (|15|) and (|22|) , constructed so as to preserve 
this asymptotic behaviour, should provide a result for the large spacelike- 
q 2 behaviour of F v (q 2 ) whose only model dependence is in the confinement 
mechanism. 

Such a model will provide an extrapolation to small spacelike-g 2 that is sensi- 
tive to the model parameters. The DSE framework provides an interpretation 
of these model parameters in terms of qualitative and quantitative features of 
the quark-quark interaction in the infrared. 

To evaluate Eq. ([!]) it is convenient to work in the Breit frame with 

g=(0,0,-g,0) and P = (0, 0, ±g, i^m\ + \q 2 ) (41) 
in which case, with k = k(sm(3 sin 9 cos 0, sin/? sin 9 sin (/>, sin j3 cos 9, cos j3), 

k ■ q= — kq sin/5 cos6> (42) 
k ■ P = \k q sin (3 cos 9 + ik\Jm 2 + \q 2 cos (3 (43) 

and one is left with three integrals to evaluate, one radial and two angular. 

The definition of confinement employed herein ensures that, for all spacelike- 
q 2 , the integrand is regular and hence that the integrals can be evaluated using 
straightforward Gaussian quadrature techniques; i.e., there are no endpoint or 
pinch singularities. 

4-1 Fitting the Parameters 

Equations (|T4"D, ( |15D and (^) provide a six parameter model [C, m, b . . . b 3 ] 
of the nonperturbative dressed-quark substructure of the pion based on DSE 
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studies. These parameters are fixed herein by requiring that the model repro- 
duce, as well as possible, the following experimental values of the dimensionless 
quantities: 

f TTl 2 

-rrr = 0.423, f v r v = 0.318 , —^ = 0.396; (44) 
{qq) 3 (qq) 3 



the dimensionless it-it scattering lengths (see Refs. pO|j3T[] for a discussion): 

(45) 



a° = 0.26 ±0.05 , = -0.028 ± 0.012 



a\ = 0.038 ± 0.002 , a° = 0.0017 ± 0.0003 ; 

and a least-squares fit to F^q 2 ) on the spacelike-g 2 domain: [0,4] GeV 2 . The 
fitting procedure was performed using the expression for the pion decay con- 
stant, f n , given in Eq. fl32]) , 



(qq)^ = (46) 

( a 2 \" ( A 2 \~ a 3 Au r v ( m 

- I In — I lim In - 9 UV —r / ds s (cr s (s) - 

*^oc \ A 2 nnn / 4tt 2 J \ W 



A qcdJ ^ 2 uv^\ a qcdJ 4tt 2 J V s + m 2 



s = p 2 and Aqcd — 0.20 GeV, following Ref. [^] but with a = 1 rather than 
the anomalous dimension, d, because the associated ln[p 2 ]-corrections have not 
been included in Eqs. ( |I~4| ) and (|T5D, and the expressions for and 
rv given in Ref. |30| . Using Eq. (|l^) in Eq. (|46|) one obtains 

fe) ^_ (2fl) l(,_g_)_LA_, (47) 

which is independent of 62 for arbitrarily small but nonzero A. 
Following this procedure one obtains 



C = 0.0422 , fn — 0.0111 , 

60 = 0.135 , 61 = 2.48 , 6 2 = 0.502 , 63 = 0.168 . (48) 

The mass scale is set by requiring equality between the percentage error in f n 
and rv, which yields D = 0.132 GeV 2 . 

The results obtained are robust in the sense that they are not sensitive to the 
functional form of the quark propagator. Different functional forms, provided 
they ensure confinement in the manner described herein, whose parameters 
are allowed to vary in order to provide a good fit to pion observables, yield a 



15 



numerical form for the quark propagator that is, pointwise, approximately the 
same on the integration domain explored by pion observables. Since the quark 
propagator implicitly constrains the gluon propagator via Eq. (||), the robust 
nature of the study makes it plausible that pion observables can be used to 
constrain the dressed-gluon propagator. 

To make the connection with the gluon propagator explicit is computationally 
more intensive, with a calculation of the pion form factor requiring a solution 
of Eq. (|j) off the real-spacelike axis. A first step in that direction is made 
in Ref. ||25|| , in which a one-parameter model dressed gluon propagator is 
employed in the calculation of all the observables considered herein, except for 
the form factor. That study illustrates that all of the parameters in the quark 
propagator are correlated via the single parameter in the gluon propagator and 
that one can fit the available pion data to the same level of agreement with only 
one parameter, in addition to the current quark mass. The quark propagator 
determined by this gluon propagator agreed well in form and magnitude with 
the parametrisation employed herein evaluated with the best-fit parameters. 



4-2 Results for F 7T (q 2 ) and Other Observables. 



In Table [l] the low-energy physical observables calculated with the parameter 
set of Eq. ( pE8D are compared with their experimental values. In this table 
m ave _ _|_ m ^ j2 anc [ m ^ j s obtained from: m 2 f 2 = — 2 m^ (qq) ^- Each 
of the calculated quantities tabulated here was evaluated at the listed value of 
m^; i.e., the chiral limit expressions for f n , Eq. (|32"D, r n , Eq. (A21) of Ref. |30|j , 
and (7,r0 77 , Eq. (|37|), were not used, however, the finite pion mass corrections 
are less than 1% for each of these. The experimental values of m ave , m n and 
g nll are extracted from Ref. ; r n from Ref. |33[ 



are discussed in Ref. 



5!| ; r n from Ref. |33[ ; and the scattering lengths 
The value of (qq) is that typically used in QCD 
sum rules analysis and the fitting error allows for deviations of 50%. The results 
presented in the table support the notion that, in an expansion of the pion 
mass in powers of the current-quark mass, the leading order term dominates. 



The results presented herein are representative of DSE studies ||22| . [25 



In Fig. H] the five n-n partial wave amplitudes associated with the scattering 
lengths in the table, calculated using the formulae in Ref. P0"| , which do not 
take final-state ir-ir interactions into account, are plotted. They are in reason- 
able agreement with the data up to x ~ 3, which corresponds to E ~ Am^. 



The form factor, i^g 2 ), at small spacelike-g 2 is shown in Fig. |3] and for larger 
spacelike-g 2 in Fig. |j. Given that the "experimental" point at q 2 = 6.3 GeV 2 , 
measured in pion electroproduction P7J, depends strongly on the model used 
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Table 1 

A comparison between the low-energy tt observables calculated using the parameters 
of Eq. (E^) and their experimental values (see text for sources) 





Calculated 


Experiment 


u 


0.0837 GeV 


0.0924 ± 0.001 


1 

-W<3ev» 


0.221 


0.220 ± 0.050 


ave 
m iGeV 2 


0.0057 


0.008 ± 0.007 


m n 


0.132 


0.138 




0.595 fm 


0.663± 0.006 




0.498 (dimensionless) 


0.500 ± 0.018 


U 


0.191 


0.26±0.05 


4 


-0.0543 


-0.028 ± 0.012 


2ag - 5ag 


0.654 


0.66 ± 0.12 


a\ 


0.0380 


0.038 ± 0.002 


4 


0.00170 


0.0017 ± 0.0003 




-0.000286 





to separate strong and electromagnetic effects, the model agrees well with the 
experimental data. 



Collectively these results indicate that the phenomenological DSE approach 
employed herein, which describes the pion as a bound state of dressed quarks 
interacting via nonperturbative gluon exchange, provides a concise, uniformly 
good description of the properties of the pion. This admits the interpretation 
that, away from resonances, the nonperturbative dressed-quark core of the 



pion is its dominant determining characteristic, as argued in Refs. [27,30 . 



5 Asymptotic behaviour 



A simple way to analyse the behaviour of F 7T (q 2 ) at large spacelike-g 2 is to 
rewrite Eq. fll]) in terms of the Bethe-Salpeter wave- function defined in Eq. 



A,(P + q, -P) 



(49) 



2N C 

J 7T 



d 4 k 



trc 



17 




Fig. 2. A comparison between the tt-tt partial wave amplitudes calculated using 
Eq. ( |48| ) (solid line) and experiment (the data are taken from Fig. 4 in Ref. |30| ). 
The dashed line is the current algebra prediction of Ref. [M]. This is identically 
zero for T®(x) and T| (x). The dimensionless variable x = E 2 /(4m 2 ) — 1 is defined 
so that threshold is at x = 0. Note that x = 3 corresponds to E = Am^. 



Asymptotic freedom, which is represented up to ln[p 2 ]-corrections by Eqs. fllOl) 
and (^), entails that the behaviour of F n (q 2 ) at large spacelike-g 2 can be 
obtained from Eq. (^) with 



«^ a,ld J * p|sii i7TAikl' (50) 

where A 1; A 2 and M c are characteristic parameters [the value of which is 
unimportant but typically [[27]] A x ~ 500 MeV ~ A 2 and M c ~ 220 MeV] 



and T fl (pi,p 2 ) ~ 7 M . It follows from Eqs. (p0|), ( f49|) and (^(J) that, at large 
spacelike-g 2 , 

F n (q 2 ) oc J Xn(k) xAk ~ \q) ■ (51) 
Making use of the approximation 
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F(q 2 ) 



1.00 



0.75 



0.50 



0.25 



0.00 




0.2 

q 2 (GeV) 2 



Fig. 3. The impulse approximation to F n (q 2 ) calculated using the parameters in 
Eq. (|48|). The experimental data are from Refs. (3(| (circles) and [53] (crosses). 



I \ d9 sin 2 9 



f(k 2 + p 2 — 2kpcos6 
k 2 + p 2 — 2kpcos9 



k 2 p l 



9(p 2 - k 2 



(52) 



(53) 



which is often used in the analysis of the asymptotic behaviour of Eq. (|]) and 
is very good for large spacelike-g 2 [j38| , one obtains 



'FL(x) 



QCXXn(x), 



(54) 



where x = q 2 /^. The solution of this equation, which satisfies the boundary 
condition F n (q 2 = oo) = 0, is 

y 

(55) 



F n (x) oc CiXv{x) + J dyx' n (y) J dzz\A z ) 
where C\ is an undetermined constant. 



One observes from Eq. ([55]) that in impulse approximation the asymptotic 
form of the elastic form factor depends on the Bethe-Salpeter amplitude of 
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q 2 F,(q 2 ) 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 




Fig. 4. The impulse approximation to q 2 F 7T (q 2 ) calculated using the parameters in 
Eq. (|48|). The experimental data are from Refs. |35| (crosses), |3(| (diamonds) and 
[37 1 (circles). 



the bound state. This result indicates the failure of the factorisation Ansatz 
in the present analysis of this exclusive process. It follows because the confine- 
ment mechanism explored herein eliminates endpoint and pinch singularities 
in Eq. (0). 

A dependence of the asymptotic fall-off of F^i^q 2 ) on the pion's Bethe-Salpeter 
amplitude is also found using the light-front formulation of relativistic quan- 
tum mechanics ]39|]. In this approach the asymptotic form of F. K {q 2 ) is only 



independent of Xw if the constituent-mass of the quark is zero, when the mass- 
shell singularity dominates the integral that arises. In such an approach, how- 
ever, a constituent-quark mass of zero is a phenomenologically untenable as- 
sumption 0, with a value of ~ 210 MeV being required to fit the available 
data. 

Using the form of x-k given in Eq. (|50|) , which is a consequence only of asymp- 
totic freedom, one finds 



FM 2 



q — »oo 
OC 



In q 2 



(56) 
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Taking into account the ln[p 2 ]-corrections to Eqs. ( [L0|) and (p3|), which arise 
because of the anomalous dimension of the propagator and Bethe-Salpeter 
amplitude in QCD, would only lead to the modification ln[g 2 ] — > ln[g 2 ] 7 , where 
|7| is 0(1), in Eq. (||). 

The numerical methods used herein to calculate F^{q 2 ) have been constructed 
so as to ensure that the result is independent of the details of the numerical 
procedure for < q 2 < 20 GeV 2 . Therefore, in order to verify the result 
of Eq. (|56"D and to estimate the spacelike-g 2 at which the asymptotic regime 
of this calculation is reached, a least squares fit of the calculated results for 
1/F^q 2 ) to 

X X 2 

a ° + ai InX +a2 InX' (57) 

with X = q 2 /(l GeV 2 ), was performed on the domain 5 < X < 20. This 
procedure yielded 

a = 39.2 , a x = - 14.7 , a 2 = 1.55 . (58) 

The fitting function and the calculated results are compared in Fig. [|, which 
confirms Eq. ([56]). This analysis suggests that, with the parameter values of 
Eq. fl48|) , which are fixed by physics at spacelike-g 2 < 4 GeV 2 , the asymptotic 
term; i.e., the ln[V]/V 2 term, only becomes dominant (provides more than 
60% of the magnitude of the form factor) for spacelike-g 2 > 10 GeV 2 . This 
result is consistent with the arguments of Ref. ||; i.e., soft, nonperturbative 
physics dominates at presently accessible spacelike-g 2 in exclusive processes. 

Equation (^) follows because of the constraints that the realisation of con- 
finement explored herein place on the fermion propagator. This mechanism is 
not an essential element of the DSE framework. The invalidity of Eq. ([55]), 
and hence Eq. fl5"6|), would not compromise the quantitative results discussed 
in the preceding sections. 



6 Summary and Conclusions 



The impulse approximation, illustrated in Fig. p], has been used to calculate 
F n (q 2 ) at spacelike-g 2 ; i.e., away from resonance contributions. The impor- 
tant elements of this calculation are: the dressed quark propagator; the pion 
Bethe-Salpeter amplitude; and the dressed quark-photon vertex, which follow 
from realistic, nonperturbative Dyson-Schwinger equation studies. This entails 
that in this calculation each of these elements behaves at large spacelike-p 2 
in the manner prescribed by the renormalisation group in QCD, up to ln[p 2 ]- 
corrections, and at small spacelike-p 2 in such a way as to ensure confinement. 
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5.0 



10.0 



15.0 



20.0 



150.0 




q 2 (GeV) 2 

Fig. 5. The upper panel compares F 1T {q 2 ) with the reciprocal of the fitting function. 
The a2 term provides more than 60% of the value of F n (q 2 ) only for q 2 > 10 GeV 2 . 
The lower panel compares 1/F n (q 2 ) (plotted points) with the fitting function of 
Eq. 



In addition, the same approach has been used to simultaneously calculate: /„-; 
m n', (<?<?); r n ; the 7r° — > 77 decay width; the 7T-7T scattering lengths: a[j, ajj, a\, 
a°, a 2 ,; and the associated partial wave amplitudes. 

The calculated results for all quantities agree well with the data, as discussed 
in Sec. [4.2| . This supports the contention that the confined, nonperturbative 
"dressed-quark core" is the dominant determining characteristic of the pion, 
away from resonance contributions. Particularly interesting in this context 
is the calculation of the 7r° — > 77 decay width, Sec. |3|. The correct chiral 
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limit value for this decay width is obtained independent of the details of 
the modelling of the quark propagator. This is an illustration of the man- 
ner in which the Abelian anomaly manifests itself in the Dyson- Schwinger 
equation approach. This calculation indicates that soft, nonperturbative con- 
tributions dominate the exclusive electromagnetic form factor for spacelike- 
q 2 < 10 GeV 2 . 

The Dyson-Schwinger equation approach provides a phenomenological frame- 
work in which to relate experimental observables to the qualitative and quan- 
titative features of the effective quark-quark interaction in the infrared, which 
is an important but presently unknown aspect of QCD. There have been at- 
tempts to calculate the gluon propagator directly: some by solving the Dyson- 
Schwinger equation for the gluon vacuum polarisation |[15HT7|; and others via 
lattice simulations |4(|, which are presently qualitatively and quantitatively 



unreliable for spacelike-g 2 < 1 GeV 2 . The present study suggests that the 
framework explored herein can be a valuable complement to such analyses. 
More precise measurements of F^{q 2 ) for spacelike-g 2 > 1 GeV 2 would be 
useful in this regard. 
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